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Signatures of photon localization are observed in a constellation of transport phenomena which 
reflect the transition from diffusive to localized waves. The dimensionless conductance, g, and the 
ratio of the typical spectral width and spacing of quasimodes, 5, are key indicators of electronic 
and classical wave localization when inelastic processes are absent. However, these can no longer 
serve as localization parameters in absorbing samples since the affect of absorption depends upon 
the length of the trajectories of partial waves traversing the sample, which are superposed to create 
the scattered field. A robust determination of localization in the presence of absorption is found, 
however, in steady-state measurements of the statistics of radiation transmitted through random 
samples. This is captured in a single parameter, the variance of the total transmission normalized 
to its ensemble average value, which is equal to the degree of intensity correlation of the transmitted 
wave, K. The intertwined effects of localization and absorption can also be disentangled in the time 
domain since all waves emerging from the sample at a fixed time delay from an exciting pulse, t, are 
suppressed equally by absorption. As a result, the relative weights of partial waves emerging from 
the sample, and hence the statistics of intensity fluctuations and correlation, and the suppression of 
propagation by weak localization are not changed by absorption, and manifest the growing impact 
of weak localization with t. 

PACS numbers: 42.25.Dd, 42.25.Bs, 05.40.-a 



OVERVIEW 



Localization 0,[ 



U, |5| is inherently a wave interfer- 
ence phenomenon If | and may therefore occur for both 
classical and quantum me- 

chanical waves. It may be considered for all manner 
of excitations including electromagnetic, sound and ul- 
trasound, acoustic and optic phonons, surface plasmons, 
polaritons, electrons and atoms. This article will discuss 
signatures of electromagnetic localization in the statis- 
tics of radiation reflected from and transmitted through 
random samples. 

The study of transport has been enriched by the many 
differences and similarities between electronic conduction 
and electromagnetic wave propagation. Differences in the 
types of wave and their scattering interactions affect the 
types of measurements that can be carried out, as well 
as the possibility and ease with which waves can be lo- 
calized. Whereas classical waves may be phase-coherent 
in large static samples since wavelengths greatly exceed 
atomic dimensions, electronic samples must be cooled to 
ultra-low temperatures in order to suppress dephasing 
even in micron-sized samples. The value of the transport 
mean free path, £, in which the direction of propaga- 
tion is randomized, is determined by the specifics of the 
wave interaction in each material system. However, the 
character of classical aird quantum propagation in meso- 
scopic systems, in which the wave is temporally coherent 
throughout the sample 0, 0, 0, , is strikingly sim- 
ilar on length scales greater than both the wavelength 
and the mean free path. For both quantum and classi- 
cal waves, the character of transport depends upon the 



closeness of the wave to the threshold of the localization 
transition separating extended and localized waves. 

loffe and Regel proposed that propagating waves can- 
not exist if scattering occurs on a length scale smaller 
than X/2tt = 1/k, in which the amplitude of the wave 
changes. This suggests that propagation ceases when 
scattering is strong enough that, £ < 2tt/X, or k£ < 1 
1^. The interference of waves that follow the trajecto- 
ries in opposite senses around a closed loop enhances the 
return of the waves and thus suppresses transport. In 
weakly scattering three-dimensional systems, the proba- 
bility that a meandering wave trajectory will loop back 
to a typical coherence volume, Vc'--"(A/2)"^,is~l/ {k£)'^ . 
The impact of wave interference on average transport is 
therefore small when k£ ^ 1. The ensemble average of 
the spread of the square amplitude of the wave is then 
diffusive and corresponds to the random walk of electrons 
or photons. However, when k£ ^ 1, the probability that 
waves will return to a coherence volume is substantial. 
Interference then impedes the escape from this volume. 
This consideration leads again to the loffe- Regel criterion 
for localization, k£ 1. 

Measurements of k£ therefore give an indication of 
the closeness to the localization threshold. Values of k£ 
can be determined from measurements of the width of 
the peak of enhanced backscattering, which is ^ l/k£ 
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This peak has a triangular cusp with 
a maximum enhancement of retro-reflection over the dif- 
fuse background of a factor of 2 in the diffusive limit 
[lil Il9,. 20], and a smaller enhancement due to recur- 
rent scattering events when the localization threshold is 
approached 12111, as seen in Fig. 1. The reflected enhance- 
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FIG. 1: Coherent backscattering of light measured, with two 
different values for the transport mean free path i. The typ- 
ical angular width varies as Xji. Narrow cone: a sample of 
BaS04 powder with £/A = 4; broad cone: TiOg sample with 
IjX = 1. The inset confirms the triangular cusp predicted by 
diffusion theory, and also shows that the maximum enhance- 
ment factor is lowered for the sample with small value of IjX 
PH . (with kind permission of the authors) 



ment peak is the Fourier transform of the point spread 
function of the incident beam [2^. 

The closeness to localization in diffusive samples can be 
ascertained from a determination of I from measurements 
of the scaling of total transmission [i^. Measurements 
of optical transmission through a wedge-shaped sample 
of random titania particles dispersed in Polystyrene are 
shown in Fig. 2 "23]. In the diffusive limit, the ensemble 
average of the spatial intensity distribution corresponds 
to the solution of the diffusion equation for the intensity. 
When absorption is absent, the average energy density 
falls linearly with depth inside the sample, corresponding 
to transmission which falls inversely with sample thick- 
ness, L. When absorption is present in a diffusive sample, 
the ensemble average of both the energy density versus 
depth and the transmitted flux versus L falls exponen- 
tially with attenuation length, La — y/Dra, for L > La, 
where D is the diffusion coefficient and is the absorp- 
tion time 23]. The variation of transmission as 1/L is 
seen in the log-log plot in Fig. 2a for L < La, while 
the exponential variation of transmission for L > La is 
exhibited in the semi-log plot in Fig. 2b. Transmission 
also falls exponentially as a result of wave localization in 
strongly scattering samples. 

In three-dimensional samples, the diffusion coefficient 
is given by, D = ^ve^, where ve is the transport veloc- 
ity li^l . The transport velocity is reduced when the wave 
is resonant with the constituents of the sample, and ex- 
hibits dips near Mie resonances with spherical dielectric 
particles [13, 113, Ell • The transport mean- free-path 
can be found once D and ve has been determined. The 
diffusion coefficient can be determined dir ectly from mea- 
surements of pulsed transmission [23, 123, y^, 1^ 1^ 



FIG. 2: Scale dependence of optical transmission through a 
wedge-shaped sample of random titania particles dispersed in 
Polystyrene. The curve gives the fit of the envelope of the 
data to diffusion theory. Transmission varies as (a) 1/L for 
L < La and as (b) exp(— L/La) for L > La = 112 ± 5 m. 
From Ref. 

or from the correlation function of the field or intensity 
with frequency shift [H El El H E^ . The field corre- 
lation function is the Fourier transform of the ensemble 
average of the response to a short incident pulse SI ■ 
The diffusion coefficient can also be determined in diffu- 
sive systems from the intensity correlation function with 
delay time in a medium in which the constituents are in 
random motion, such as in colloidal samples [s^lssj . The 
temporal correlation function reflects the average path 
length distribution in the medium, which is proportional 
to the average temporal evolution of transmitted pulses. 

For bounded samples, in the absence of inelastic pro- 
cesses, however, any marker of localization must also de- 
pend upon the sample geometry, in addition to i and A, 
since reflection from the sample boundaries also returns 
the wave to points within the sample. When inelastic 
scattering is present, either the impact of inelastic scat- 
tering must be explicitly accounted for or a propagation 
parameter must be found, which still captures an essen- 
tial aspect of localization and thereby directly expresses 
the closeness of the wave to the localization threshold. 

Spectra of transmitted intensity for extended optical 
E3 and microwave [s^ E^l radiation, such as the mi- 
crowave measurements shown in Fig. 3a ;3Qj ■ are com- 
posed of many overlapping lines which result in spectra 
without isolated peaks. The wave inside the medium, 
as well as each of the quasimodes, is extended in space. 
Within the frequency range in which the wave is local- 
ized, however, distinct resonances are observed in the 
transmission spectrum [3, 0|i as seen Fig. 3b. This 
spectrum of microwave transmission of localized waves 
was measured in a sample of alumina spheres at low 
density, contained in a copper tube, which was cooled 
to the temperature of liquid nitrogen to reduce absorp- 
tion losses in reflection from the tube. When the wave 
resonantly excites the sample, its amplitude may be ex- 
ponentially peaked within the sample [i^ E3i Ej| • The 
relatively high intensity in the interior of the sample, as 
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FIG. 3: Isolated peaks in (b), as opposed to (a), are an indica- 
tion of photon localization. The variance of relative intensity 
is high in spectra with isolated peaks. From Ref. [39l |. 



compared to near the boundaries, leads to long decay 
times for the excited mode and, consequently, to distinct 
narrow spectral lines. Transmission via localized quasi- 
modes is particularly high when intensity is peaked near 
the middle of the sample 4^ ,43j| . Transmission is also 
high when two or more quasimodes occasionally overlap 
spectrally leading to a multi-peaked intensity distribu- 
tion exten ding through the sample, known as necklace 
states 0, 0rE3 ■ Because these short-lived states are 
spectrally broad, they dominate average transmission in 
ID system. Off resonance, the overall energy density falls 
exponentially within the sample. 

Localizing photons has proven to be a considerable 
challenge. Photon states bound to subwavelength-scale 
particles do not exist as they do for electrons bound to 
ionic cores. In addition, s-wave scattering does not exist 
for photons, as it does for electrons. The electromagnetic 
scattering strength due to p-wave scattering vanishes in 
the long- wavelength, Rayleigh limit, as a^/A^, where a 
is the particle radius. As a result, it is not possible to 
enhance scattering by packing many small scattering par- 
ticles into volumes comparable to the wavelength. Scat- 
tering for spherical particles is peaked at Mie resonances. 



when the particle radius is comparable to A. But then it 
is difficult to create strong scattering with £ as short as A, 
since the maximum strength of scattering from individual 
particles is achieved only in dilute systems when scatter- 
ers are separated from one another. Because photons are 
uncharged, however, a pure Anderson localization tran- 
sition between diffusing and localized waves might occur 
without the complication of the Coulomb interaction. If 
this could be observed, it would be possible to deter- 
mine the critical exponents of the Anderson transition 
[iil li^ . In addition, the statistics of classical trans- 
port can be determined in large ensembles of equivalent 
samples. To date, photons have only been localized in 
three-dimensional systems with residual long-ran ge o rder 
m 113, in one- [El 112, two- js^Land quasi-lD 113 
samples, and in layered samples [So. I57ll58l| . 

John pointed out that waves can be localized in 
nearly periodic saniples in the frequency range of the 
photonic band gap 1^ that forms in the associated 
periodic structure when the index of modulation is suf- 
ficiently large. Though the density of states vanishes in 
large periodic structures, the introduction of disorder cre- 
ates localized states in the spectral region of the photonic 
band gap. The low density of states within the spectral 
range of the photonic band gap and the narrow, spec- 
trally isolated lines associated with the long residence 
time of light on resonance with these states lead to spa- 
tially localized modes. Measurements of large microwave 
intensity fluctuations within the low frequency band gap 
of a cubic wire structure have indicated that localiza- 
tion occurs for intermediate filling fractions of metallic 
spheres floated within the structure 0, js^l . 

The relationship between average transmission and 
fluctuations in transmission lies at the core of the scaling 
theory of localization 0, 0,E3- 1^, it has been hy- 
pothesized that the relationship can be expressed via a 
single parameter 0|, which relates fluctuations in trans- 
mission, T, to the ensemble average of the logarithm of 
transmission, which is normally distributed. The vari- 
ance of the logarithm of the transmission in samples of 
length, L, was conjectured to be, = j/L, where 7 is 
the Lyapunov exponent, which is the inverse of the local- 
ization length, 7 = 1/^ = — InT/i, for localized waves 
0. Recently strong deviations from the single parame- 
ter scaling hypothesis have been predicted in the tail of a 
band in which the density of states is low [HsllH^ IgsL l66| ■ 
In addition to the localization length, ^, the scaling of 
transmission then depend s up on a second length scale, 
£s, whenever ^ < £s |E1> EJ- This distance is the av- 
erage spacing between localized states with frequencies 
between the given frequency and the edge of the band at 
which the density of states vanishes. The breakdown of 
single parameter scaling has also been demonstrated in 
2D samples [slElES- 

There has been considerable interest in propagation 
and localization in quasi- ID samples. These samples 
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are locally three-dimensional, but have reflecting trans- 
verse boundaries with length, L, considerably exceeding 
the sample's transverse dimensions. Since all modes are 
mixed in such systems, the statistics of propagation are 
the same at ai iy p oint on the output surface. Random 
matrix theory [7ll is well suited to such systems and 
has provided nonperturbative results that supplement 
the findings of diagrammatic perturbation theory in the 
weak scattering limit 0,E3|- the same time, measure- 
ments of field, intensity, and total transmission in the 
frequency and time domains can be performed on ran- 
dom ensembles of quasi- ID samples. Because of reflection 
at the transverse boundaries, the wave is inhibited from 
straying far from a point in the transverse directions. 
The number of returns to a typical coherence volume, 
Vc ~ (A/2)'^, increases superlinearly, as L increases. As 
a result, the wave can always be localized as the sample 
length is increased while the cross sectional area remains 
fixed. This makes it possible to obtain a comprehensive 
statistical picture of propagation and localization in sam- 
ples with well-defined statistical parameters in diffusive 
and localized samples [t^. 

Key properties that lend themselves to a universal de- 
scription are the average values of transmission and con- 
ductance and their scaling, and the correlation and vari- 
ance of local and spatially averaged fiux or current. A 
variety of parameters capture the changing character of 
diverse aspects of wave propagation in the localization 
transition. The relationship between these phenomena 
and the parameters that characterize them are explicit in 
the absence of inelastic processes, since transport may be 
described in terms of a single parameter near the center 
of a band. Indications of a breakdown of single parame- 
ter scaling in regions in which the density of states is low 
have not been found in measurements in locally 3D sam- 
ples in the slab or quasi-lD geometries. When absorp- 
tion is present, however, scaling cannot be described in 
terms of a single parameter, and the resultant exponen- 
tial reduction in transmission is not related to increasing 
localization of the wave. Though the conductance can 
then no longer serve as a localization parameter, mea- 
sures of correlation and fluctuations of transmitted flux 
may still be expressed through a single parameter which 
is a reliable indicator of the nature of wave transport. 

We will find that a clearer picture of wave propagation 
can be achieved once the effects of absorption and weak 
localization are disentangled. This can be done by trans- 
forming from the frequency to the time domain, since 
absorption has the same impact on all paths at a given 
delay from an exciting pulse. Since the distribution of 
paths is not changed by absorption, relative fluctuations 
are unchanged at a flxed delay time t from an exciting 
pulse. The temporal evolution of the wave can also be 
explored by studying its Fourier transform, the fleld cor- 
relation function with frequency shift. 

In the absence of inelastic processes, the suppression 



of transport due to the enhanced return of the wave to a 
coherence volume may be described in terms of the av- 
erage value of the dimensionless conductance g, which 
is the scaling parameter 0, Q of the Anderson local- 
ization transition, g = G/ (e^/i). In a particular sam- 
ple realization, the dimensionless conductance equals the 
transmittance T, which is the sum of transmission coeffi- 
cients over all N input transverse modes, a, and outgoing 
modes, b, T ^ ^abTab |llEi,E3- The conductance is 
naturally measured in electronic experiments since the 
relative phase of the incident transverse current modes 
fluctuates randomly over the measurement time to yield 
the incoherent sum of transmission coefficients. In con- 
trast, studies of classical waves are often carried out with 
coherent sources, so that it is natural to measure either 
the field or fiux transmission coefficient for a single inci- 
dent mode to a single outgoing mode or the total trans- 
mission which is the transmission coefficient for the net 
transmitted flux for a single incident transverse mode. 

In quasi-lD samples, the number of transverse modes 
in which two transverse wave polarizations are distin- 
guished is = 27rA/A^, where A is the sample cross 
section. In the diffusive limit, g ^ 1, the dimensionless 
conductance can be approximated by g — N£/L. This 
is equivalent to Ohm's law for diffusing particles, since 
g = (2Tr£/X^) {A/L). To order 1/g, the conductance is re- 
duced from the expression above by a factor of (1 — 1 /3g) 
[TSf . Scaling of average transport is given in terms of 
the single parameter, g, and the localization threshold is 
reached at g ~ 1 , at which point g is substantially re- 
duced by renormalization due to weak localization. Since 
the value of g falls at least as fast as it does for diffusing 
particles, it is possible to localize waves in static systems 
by choosing sample lengths greater than the localization 
length, ^ = N£. This follows the analysis by Thouless 
that electrons can be localized in sufficiently long wires 
when the temperature is low enough that dephasing can 
be neglected Q]. 

The relationship of g to localization may be seen from 
its connection to the Thouless number, S 0, which 
may be expressed as the ratio of the typical width and 
spacing of resonances of a random ensemble of samples, 
S — Si^/Av. The level width is the average leakage rate 
of energy from the electromagnetic quasimodes of the 
sample and is essentially the inverse of the width of the 
distribution of photon transit times through the sample 
when inelastic processes are absent. This corresponds to 
the width of the field correlation function with frequency 
shift jz^ and is given hy 5 k, D/L'^ 0. The level spacing 
is the inverse of the density of states of the sample as a 
whole, and may be written as lS.v = \/n{i')AL^ where, 
n(^) is the number of states per unit volume per unit 
frequency, and AL is the sample volume. The Thou- 
less number is a measure of localization since transport 
is then suppressed by the poor overlap of quasimodes in 
different sectors of a sample when, 5v < Aj/, and hence. 
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6 < 1. When S > 1, on the other hand, quasimodes 
overlap spectrally and transport is relatively uninhib- 
ited. Using the Einstein relation for the conductivity, 
a = {e^ / h) Dn{v) , the conductance can be written as 
G = ctA/L = (e^/Zi) Dn{v)A/L. We can then write, 
5 = 5v//^v = {D/L^) I \\ln{y)AV) = Dn{v)AIL = 

Gl (eV/i)=ff 0. 

The degree of correlation can be identified with 1/(5, 
since, for monochromatic excitation, the number of quasi- 
modes which are substantially excited is approximately 
the number of modes within the mode linewidth, which 
is 5v / IS.V = 5. Since the state of the wave in the sample 
is approximately specified by the amplitudes of these 5 
modes, it may be specified approximately by 5 parame- 
ters. As a result, the degree of correlation of intensity is 
approximately, 1/(5, or equivalently, 1/g. 

In the diffusive limit, 5 3> 1, far from the localiza- 
tion threshold, the change in average conductance and 
the degree of long-range correlation in the flux due to 
wave interference are both of order 1/g. Nonetheless, 
fluctuations in spatially averaged flux are dramatically 
enhanced over the values that would be obtained if long- 
range correlation were absent and the flux in each of 
the N transmission modes could therefore be assumed 
to be statistically independent. If this were the case, 
the variance of total transmission normalized to its en- 
semble average value, Sa = Ta/{Ta), would be given by 
var(sa) = ^/N, while the variance of the conductance 
normalized to its ensemble average, s — T/g, would be 
given by var(s) = l/N"^, when absorption is absent. In- 
stead, the variance of normalized total transmission over 
an ensemble is enhanced by a factor of L/£ , to 
var(sa) « {l/N){L/i) « 1/ff HI US El S El 
The variance of the normalized dimensionless conduc- 
tance is enhanced by a factor of (Ljt)^ , to give var(s) ~ 
(1/7V2) {Litf w l/g2 [ll[llE3- This leads to a fixed 
value for the variance of the dimensionless conductance, 
var(T) = (7^var(s) « 1, where, T = s x g. Calculations, 
in accord with measurements, give var(T) — 2/15, in- 
dependent of the size or local scattering strength of the 
sample 0,0,1^. This is termed universal conductance 
fluctuations. 

The neglect of absorption in calculations of propaga- 
tion makes it possible to highlight specific properties re- 
lated to weak localization and mesoscopic fluctuations. 
However, this is only rarely justified in practice. Indeed, 
the attempt to deal with the reality of absorption has 
sharpened our understanding of the nature of localiza- 
tion and has provided clear signatures of localization. 



ENSEMBLE AVERAGE TRANSMISSION AND 
REFLECTION 

Early discussions of photon localization were stimu- 
lated by the analogy between electron and photon trans- 
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port. The impact of localization on the scaling of trans- 
mission was represented as a renormalization of local 
scattering parameters such as the mean free path or dif- 
fusion coefficient js^. It was conjectured that the renor- 
malization of optical transport parameters would be cut 
off by absorption jus t as electron renormalization is cut 
off by dephasing [lOj. Microwave measurements of the 
variation of transmission with thickness in strongly scat- 
tering mixtures of metal and dielectric spheres showed 
departures from diffusive behavior that were interpreted 
in terms of a scale dependent diffusion coefficient whose 
value saturated in samples Icmger than the attenuation 
length due to absorption jsl E3 ■ 

On length scales shorter than the absorption length, 
the effective diffusion coefficient at the mobility edge is 
predicted to scale as 1/L and the transmission as 
while the transmission is expected to fall exponentially 
for localized waves 0, 0|. Such scaling was observed 
in microwave measurements [sH. E3 ) but absorption also 
leads to a more rapid fall-off of transmission in diffusive 
media j23|. Subsequently, similar scaling was observed 
in optical measurements in very weakly absorbing and 
strongly scattering samples of GaAs powders together 
with a broadening of the coherent backscattering peak 
E^. But the role of absorption might still be signifi- 
cant enough to affect these measurements in view of the 
lengthened dwell time for waves in strongly scattering 
samples E3- The closeness to the localization threshold 
in such systems can be ascertained, however, in a variety 
of statistical properties of the scattered wave, which will 
be discussed below. 

The interaction of absorption and localization was 
probed in computer simulations of the temporal varia- 
tion of the spatial distribution of acoustic energy in an 
absorbing 2D random sample that is isolated from its 
surroundings. The vibrational amplitude showed an ex- 
ponential fall-off with displacement from the single line 
which was initially vibrating, that persisted in time |90j . 
This demonstrated that though net energy transport is 
reduced by absorption, the wave is localized and does 
not spread throughout the sample. This is understand- 
able since partial waves following a loop in two opposite 
senses are in phase, when these return to the initial scat- 
tering center. Thus the impact of weak localization is 
undiminished for the surviving partial waves. 

The time dependence of flux within a medium can be 
treated separately from the loss of energy due to dis- 
persed absorption. The statistical impact of absorption 
at any time can be accounted for by simply multiplying 
all waves by the factor exp(— t/ro). This was demon- 
strated in measurements of field spectra of microwave 
radiation transmitted through a random ensemble of ab- 
sorbing dielectric quasi- ID samples 0. The measured 
probability distribution of intensity was shown to be in 
accord with a single-parameter distribution, which is the 
measured variance of the intensity, as will be discussed 
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below. It was further shown that the impact of absorp- 
tion could be undone in a statistical sense by simply 
restoring the magnitude of the amplitude of all partial 
waves emerging from the sample to the value these would 
have had if absorption were absent without altering the 
wave. This is done by transforming field spectra into the 
time domain by multiplying each spectrum by a Gaussian 
spectrum associated with the Gaussian temporal profile 
of the incident field amplitude. This product spectrum is 
then Fourier transformed to give the temporal response 
to an incident pulse with the prescribed bandwidth. The 
field in the time domain is then multiplied by exp {t/2Ta) 
and transformed back to the frequency domain. The field 
transmission coefficient is then obtained by dividing by 
the Gaussian field spectrum. The complex square of this 
spectrum is then computed to give the intensity spec- 
trum from which the probability distribution of trans- 
mitted intensity is determined. The resultant probabil- 
ity distribution of steady-state intensity was found to be 
in excellent agreement with the prediction for samples 
without absorption [H El l^l IS 1^3 • The distribution 
is then given in terms of the single parameter, g. This 
demonstrates that statistically the impact of absorption 
at any time is simply to reduce the wave amplitude. Ab- 
sorption thus has a straightforward affect on transport in 
the time domain, and does not affect weak localization 
or the associated renormalization of scattering parame- 
ters at a specific delay time t. This is confirmed by our 
finding that samples with identical scattering strengths 
but different values of absorption have the same degree 
of correlation at a given delay time. Indeed, correlation 
is not suppressed by the increasing role of absorption for 
longer time delays, but is rather enhanced with increas- 
ing time delay. The reduction of steady-state intensity 
correlation in the presence of absorption is thus seen to 
be a result of the decreasing weight of longer paths with 
increasing absorption. Finally, this analysis shows that 
the diffusion coefficient is independent of the value of the 
absorption coefficient. 

The role of weak localization increases with time since 
the probability that a wave trajectory will cross itself 
increases. We expect therefore that transport will be 
increasingly suppressed with time even in nominally dif- 
fusive samples with g > 1. An alternate picture of dy- 
namic propagation would be that the initial pulse ex- 
cites a distribution of quasimodes in the sample. The 
decay rate of transmission would then decrease in time 
since the amplitude of longer-lived quasimodes, which 
are more strongly peaked near the center of the sam- 
ple, would increase relative to that of shorter-lived res- 
onances. Averaged over an ensemble of random sample 
realizations, the distribution of leakage rates for quasi- 
modes would be a continuum and would give rise to a 
continuous reduction of the leakage rate from the sam- 
ple. Both of these pictures of pulse dynamics are at vari- 
ance with the particle diffusion picture that gives the 



ensemble average of intensity as a discrete sum over dif- 
fusion modes. The decay rates of the diffusion modes are 
given by 1/t„ = n^n^D/ (L -f 2zo)'^, where n is a posi- 
tive integer and zq is the boundary extrapolation length 
0) HE US ■ After a time ri , the intensity distribution 
would settle into the lowest diffusion mode and would 
decay at the rate, 1/ri — ■k'^D/ [L -f 2zo)^. 

Non-exponential decay has been observed in acous- 
tic scatter ing in reverberant rooms [s^i and solid 
blocks 99, 100] as well as in microwave scatter ing in cav - 
ities whose underlying ray dynamics is chaotic 10llll02l| . 
The decay rate of electronic conductance has been pre- 
dicted to fall as a result of the increasing weight of 
long-lived, narrow linewidth states |l0.4 Il04j . The lead- 
ing correction to the diffusion prediction for the elec- 
tro n sur vival probability Ps{t) was calculat ed b y Mir- 
lin |l05j using the supersymmetry approach |l06| to be, 

-lnPs(i) = (</Ti)(l-V2^Vi)- 

In pulsed microwave measurements, in quasi- ID ran- 
dom dielectric media composed of random mixtures of 
low-density alumina spheres, we find a breakdown of the 
diffusion model in the non-exponential decay of trans- 
mission j83l | . The decay rate for the transmitted pulses 
shown in Fig. 4a for different sample ensembles is seen 
in Fig. 4b to decreases in time at a nearly constant rate. 
A linear fall of the decay rate would be associated with 
a Gaussian distribution of rates of decay for quasimodes. 
We find a slightly more rapid decrease of the decay rate 
which is associated with a slower than Gaussian fall-off 
of the decay time distribution. The temporal variation of 
transmission can also be described in terms of the grow- 
ing the impact of weak localization on the dynamic be- 
havior of waves, which can be expressed via the renormal- 
ization of a time -dependen t diffusion constant or mean 
free path [H IH [lOl [Tol [109] . 

Transmission for localized waves is via resonances with 
isolated localized states or via spectrally and spatially 
overlapping states ^ Ql E^- Steady-state mea- 
surements of transmission through localized states have 
been made in stacks of glass slides [571 and in alternat- 
ing solid and porous silica layers [5g. Both overlapping 
and isolated spectral lines are observed. The presence 
of overlapping lines in the transmission spectrum of a 
layered material may have a number of causes. When 
a plane wave enters a layered sample in which isotropic 
layers are homogeneous, transmission can be described 
by a ID model. When ^ < L, then, (5 < 1, and spectral 
lines are generally distinct. However, occasional overlap 
may occur. This has a strong impact on average trans- 
mission since the spatial distribution within the sample is 
quasi-extended and the transmission is consequently high 
|46l |47| . In samples with transverse inhomogeneity, it is 
necessary to go beyond a ID model. Spectral overlap is 
then produced by broadening associated with the widths 
of the distributions of the wave vectors in the incident 
beam and of the structure factor of the layered medium 
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FIG. 4: (a) Average pulsed transmitted intensity in alumina 
samples of L = 61 (A), 90 (B and D), and 183 cm (C); Sample 
D is the same as Sample B except for the increased absorption 
by titanium foil inserted along the length of the sample tube, 
(b) Temporal derivative of the logarithm of the intensity gives 
the rate 7 of the intensity decay due to both leakage out of 
the sample and absorption. The dashed curves are the decay 
rates given by the diffusion model. From Ref. |33| . 



with transverse inhomogeneity. These additional factors 
give contributions to the fc-vector distributions which are 
separately proportional to the inverses of the width of the 
incident beam and the transverse coherence length of the 
sample. The overlap of modes is clearly seen in the com- 
plex temporal evolution of the pulse in layered porous 
silica samples 58]. Microwave measurements on a slot- 
ted single-mode waveguides packed with random slabs 
show the impact of overlapping quasimodes within the 
sample [s^ . 

Lasing is facilitated in localized modes because the ex- 
citation is able to penetrate deep into the sample with 
greatly amplified intensity when the pump laser is res- 
onant with an exponentially localized mode |57l |. Light 
is subsequently emitted into long-lived modes overlap- 
ping the region of excitation. This enhances the oppor- 
tunity for stimulated emission and therefore lower s th e 
laser threshold. In contrast, in diffusive media |llCl| . 
excitation is restricted to the surface region by multi- 
ple scattering and the lasing threshold is not substan- 
tially reduced below the threshold for amplified stimu- 



lated emission in homogeneous media Localized 
modes and low threshold lasing were observed in stacks 
of glass slides with interspersed dye solution or dye-doped 
plastic sheets ^ . Sharp- line ra ndom lasi ng spectra are 
also observed in diffusive media 

EH mi, 

in which las- 
ing is initi ated in a singl e or small number of long-lived 
modes 0; IllSl llld j . Lasing also occurs in localized 
states in nearly periodic anisotropic ID sample. This is 
observed in disordered cholesteric liquid crystals. In the 
ordered structures, the director pointed along the aver- 
age molecular axis in each layer rotates with well-defined 
pitch. A distinctive progression of band edge modes is 
then observed. But in liquid crystals in which periodicity 
has been disrupted to a lesser or greater extent, locahzed 
states form at th e band edge. Lasing then occurs in these 
states 117llll8l |. Lasing is also observed when a localized 
state is introduced into the band gap by a discontinuous 
jump in the angle of the director |ll9lll2(]| . 



STEADY-STATE STATISTICS 

In the absence of absorption, there is an explicit con- 
nection between average transport and mesoscopic fluc- 
tuations and correlation in quasi-lD samples since these 
can all be expressed in terms of a single parameter, g. 
Since this direct link between the averaged transmission 
and the statistics of transmission is lost in absorbing sam- 
ples, it is worthwhile to consider the statistics of relative 
fluctuations. We consider the probability distributions 
of transmitted intensity and total transmission normal- 
ized by their ensemble averages, Sab = Tab/ (Tab) and 
Sa = Ta/ {Ta), respectively. The relationship between the 
probability distributions of intensity and total transmis- 
sion was found by Kogan and Kaveh from random matrix 
theory jgj, 



P{Sab) = / -P{Sa)exp{-Sab/Sa) ■ (1) 

Jo Sa 

This is equivalent to the relationship between the mo- 
ments of the distributions of Sa and Sab, 



w:b)^n\{.s:), (2) 

for the n-th moment '9?| . From the relationship between 
the second moments, obtained from Eq. (2) for n — 2, 
we find, var(sab) — 2var(s£i) + 1. In the limit of weak 
scattering and vanishing absorption, the distribution of 
total transmission is found from diagrammatic and from 
random matrix theory to be |9lll92| . 



/zoo 7 
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where, 



= g In ( Vl + P/9 + VP/g 



(3) 



(4) 
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This leads to a connection between fluctuations and aver- 
age transmission, var(sa) = 2/85. Since the probabihty 
distributions of intensity and total transmission are given 
in terms of the single parameter, in this limit, these 
distributions may equally be expressed as functions only 
of var(sa). 

The distributions of Sa and Sab in three quasi-lD 
samples composed of randomly positioned Polystyrene 
spheres with different diameters and lengths in the fre- 
quency range 16.8-17.8 GHz, in which the variation in the 
mean free path is small, are shown in Fig. 5 [93lll2l| . The 
distributions broaden, and the deviation from a Gaussian 
become more pronounced, as either the sample length 
increases or the cross-sectional area decreases. The re- 
sults shown in Fig. 5 are for samples, a, 6, and c, which 
would have values of the dimensionless conductance in 
the absence of localization corrections and absorption, 
g = Ne/L, of 15.0, 9.0, and 2.25, respectively The 
straight line in Fig. 5b is a plot of exp(— Sah), which is 
the prediction of Rayleigh for the intensity distribution 
for one polarization component of the field in a sample in 
which the number of transverse modes, N, diverges and 
in which long-range correlation vanishes. The distribu- 
tion of total transmission in the diffusive limit is Gaussian 
with variance 2/3g. This is in contrast to an estimated 
value of 1/N, if long-range correlation were neglected. 

The theoretical expressions for the full distribution 
functions in Eqs. (1) and (3) are given as functions of 
g for nonabsorbing samples. We now consider the full 
transmission distributions in the presence of absorption. 
Since the decreased value of g due to absorption does not 
lead to enhanced values of intensity correlation but rather 
the reverse, it is not possible to utilize the value of g in 
Eqs. (1) and (3). But since var(sa) = 2/3(7, P{sa) may 
be expressed as a function of var(sa). When var(so) is 
calculated for the three samples in Fig. 5, and 2/3var(sa) 
is substituted for g in Eq. (4), the resulting curves plotted 
in the figure are in excellent agreement with the measured 
distributions. The measured values of 2/3var(sa) may be 
compared to the values of g calculated without absorp- 
tion which were given above: 10.2 for sample b instead 
of 9, and 3.06 for sample c instead of 2.25. The mea- 
sured values of 2/3var(sa) are larger when absorption is 
present, because absorption reduces the weight of longer 
paths. The distribution of Eq. (4 ) gi ves the exponential 
tail P(sa) ~ exp[-2sa/3var(sa)] liiliil. These results 
are in excellent agreement with measurements |93j. 

We find that var(sa) does not drop sharply with in- 
creasing absorption as does the average transmission, but 
is only weakly suppressed jij. This indicates that the 
broad distribution of intensity and of total transmission, 
which is characteristic of wave propagation in random 
systems, is preserved in absorbing media. It is of interest 
therefore to examine the scaling of the statistics of prop- 
agation. The scaling of var(sa) determined at a number 
of frequencies is shown in Fig. 6 Var(sa) increases 




FIG. 5: Distribution functions of (a) the normalized total 
transmission and (b) transmitted intensity, P{sa) and P{sab), 
respectively, for three Polystyrene samples with dimensions: 
a) d—7.5 cm, L — 66.7 cm; 6) d—5.Q cm, L — 50.0 cm; c) 
d=5.0 cm, L = 200.0 cm. Solid lines represent theoretical 
results obtained from Eqs. (1) and (3), with measured values 
of 2/3var(sa) substituted for g. The dashed line in (b) is a 
semi-logarithmic plot of the Ray leigh distribution, P{sab) = 
exp(-Sab). FromRefs. |9lll2l|. 



exponentially once it becomes of the order of unity. This 
is expected for a localization parameter, since, in the ab- 
sence of absorption, var(sa) ^ 1, when 5^1. The avail- 
ability of a measurable localization parameter makes it 
possible to definitively index the nature of propagation 
via a parameter marking the closeness to the localization 
threshold. The measurement of var(sa) makes it possi- 
ble to probe the array of factors that foster localization, 
including size, concentration, and structural correlation. 



WINDOW OF LOCALIZATION 

In the absence of absorption, localization can always 
be brought about by increasing the length of a sample, 
while holdings its cross-sectional area fixed . To induce 
localization in as short a sample as possible, strongly 
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FIG. 6: Scaling of var(sa) in samples of alumina spheres. 
Above a value of the order of unity, var(sa) increases expo- 
nentially. From Ref. ,41,). 



scattering samples should be used. Since the strength 
of scattering is peaked near Mie resonances, it is useful 
to carry out measurements in samples containing spher- 
ical scatterers at wavelengths comparable to the diam- 
eter of the spheres. In order for the resonances not to 
be washed out by correlated scattering among neighbor- 
ing spheres, samples with low density of spheres were 
studied. Microwave measurements were made in quasi- 
ID samples o f alu mina spheres at a volume fraction of 

/= 0.068 mini. The alumina spheres have a diameter 

of 0.95 cm and an index of refraction of n = 3.14. These 
are placed at the centers of 1.9-cm-diameter Styrofoam 
spheres. The sample is contained in a 7.3-cm-diameter 
copper tube with sample lengths up to L = 90 cm. An en- 
semble of 5,000 random sample configurations is created 
by rotating the tube between subsequent field spectrum 
measurements. For each sample configuration, the field 
spectrum yields the corresponding frequency variation of 
the transmitted intensity lab- The average intensity {lab) 
for L = 80 cm samples, shown in Fig. 7a, exhibits dis- 
tinct drops near Mie resonances. The resonant character 
of scattering is further indicated by the sharp peaks in 
the average photon transit time, (t), seen in Fig. 7b. 
The transit time or single-channel delay time is given by 
Tab = Sabd(j)ab/duj, whcrc (pab IS the phasc accumulated by 
the field as it propagates through the sample, from inci- 
dent chan nel a t o ou tgoing channel b, and lo is the angular 
frequency |l23l Il24l | . Since low transmission can be due 
to absorption and long dwell time can be associated with 
microstructure resonances |3 US j the average trans- 
mission and dwell time do not provide definitive measures 
of the closeness to the localization threshold. This can 
be obtained, however, from the measurement of var(scib), 
shown in Fig. 7c. For diffusive waves obeying Rayleigh 
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FIG. 7: (a) Average transmitted intensity, (J), (b) aver- 
age photon transit time, (r), (c) variance of normalized 
transmitted intensity, var(sa6), and (d) dimensionless ratio, 
{(/>'/) = /{<j}'){I), versus frequency in a quasi-lD alumina 
sample with L = 80 cm and alumina volume filling fraction 
/ — 0.068. The dashed lines i n (c ) and (d) indicate the local- 
ization threshold. From Refs. 
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statistics, va.r{sab) = 1- Lower values of var(sab), seen 
below 8.5 GHz, indicate a significant ballistic component 
in the transmitted field, whereas higher values indicate 
the presence of substantial long-range correlation. The 
horizontal dashed line in Fig. 7c represents the local- 
ization threshold, var(sa6) = 7/3, which corresponds to 
var(sa) — 2/3. In alumina samples with L = 80 cm, the 
localization threshold is crossed in a narrow frequency 
range above the first Mie resonance. Large fluctuations 
of intensity reflected in measurements of var(sab) above 
the localization threshold, seen in Fig. 7c, are a result of 
the greatly enhanced tail of the intensity distribution for 
localized radiation associated with the strong difference 
in transmission when a sample is on or off resonance. 
Within the localization window, var(saf,) > 7/3, the 
probability distribution for intensity becomes extraordi- 
narily broad 0- 

The strength of scattering cannot be substantially in- 
creased by increasing the concentration of randomly po- 
sitioned alumina spheres. When the density of alu- 
mina spheres is increased individual sphere resonances 
are washed out. This is shown in Fig. 8a in measure- 
ments of -D as a function of volume filling fraction, /, 
of nearly spherical alumina particles with approximate 
diameter 0.95 cm and index of refraction n = 2.95 [j^ . 
Measurements of the mean free path for / = 0.30 are 
shown in Fig. 8b |l2,'il |. These results yield the transport 
velocity shown in Fig. 8c, which is substantially reduced 
over a broad spectral range in agreement with CPA cal- 
culations ,126] . 



DELAY TIME STATISTICS 

In addition to the statistics of steady-state transmis- 
sion, there exists a rich statistics of dynamics, which 
is st r ongly corr elated to static aspects of propagation 
|l03L Il27l Il28l |. For a spectrally narrow pulse, the 
single-channel delay time is found to be the derivative 
of the phase of the transmitted field with frequency, 
Tab = Sabdcjiab/dio = (f)' ,123. 124] . The single-channel 
delay time can similarly be defined for the reflected field. 
The interplay between the single-channel delay time and 
transmission coefficient is expressed in the joint proba- 
bility distribution P{(j>',I). The phase is discontinuous 
with a jump of tt radians, as the null in intensity in the 
speckle pattern is traversed |l29l] . giving rise to strong 
fluctuations in at low intensities. For a given value 
of intensity / the delay time 0' is normally distributed 
with normalized standard deviation, A0'/(0') = y^Q/2I, 
whereas the second-order cumulant correlator vanishes, 
((/)'/) - {(/)') (1) = 0. For, ^ < L < ^, in the diffu- 
sive regime, in the absence of absorption, the probabil- 
ity distribution of the delay times is P(0' = cj)' / {(j)')) = 
(Q/2)\Q + {(f)' — 1)^]"'^/^, where Q « 1 in transmission 
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FIG. 8: (a) Frequency dependence of the diffusion coeffi- 
cient for mixtures of industrial-grade 0.95-cm-diameter alu- 
mina spheres within hollow polypropylene spheres for various 
alumina filling fractions, /; (b) Transport mean free path for 
sample with / = 0.30; (c) Transport velocity for f = 0.30 
obtained from (a) and (b), ve = 2,D/l. From Refs. p7i..l25,l. 



Absorption can be accounted for simply by a change in 

y^his agrees well with microwave |130t il31j and optical 
measurements. 

The statistics of dynamics of localized waves differ fun- 
damentally from those of diffusive waves since the trans- 
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mission spectrum includes narrow spikes, since S < 1. 
Long delay times for localized waves occur at peaks in 
transmission, associated with resonant tunneling through 
localized states, so that (0'/) — {4>'){I) > 0- The corre- 
lator between </)' and / thus provides a dynamical test 
of photon localization. The frequency variation of the 
ratio {4>'I) — {(j)' I) / {(j>'){I) in a quasi-lD sample of alu- 
mina spheres with L — 80 cm is plotted in Fig. 7d 
p,34j. It is unity in the diffusive limit and rises above 
unity as localization is approached. and var(sa) 

both rise appreciably only within the localization win- 
dow in Fig. 7. The localization threshold, at which 
var (sah) = 7/3, corresponds to the condition {(f>'I) =1.1, 
shown as the dotted line. The conditional probability 
distribution P/ ((/>') given / for localized waves in trans- 
mission is found to exhibit an exponential fall-off, with an 
asymmetry in the distribution, which increases with de- 
creasing /, and to hav e a n ormalized standard deviation 
of A<?!)'/(0') cx [133. The distribution, of 

the sin gle-ch annel delay time is given by random-matrix 
theory |l32j . In the localized regime, it is found to hav e 
a universal quadratic tail, P{(j)') oc It/i'P^ |l35lll36lll37l| . 
In samples of the finite length L, this algebraic tail even- 
tually crosses over to a log-normal tail, at exponentially 
long delay ti mes, Id^'l > r.^ exp(_L /f ), where Ts is the scat- 
tering time \w4 . llOSl Il38l Il39l | . There is an interesting 
dynamic signature of localization in reflection: the distri- 
bution P{(t>') differs in the two cases, a = a' and a ^ a' . 
The maximal value of P{4>') in the case a = a ' is g reater 
than that for a a' by a factor close to |l32l |. The 
effect appears only in the localized regime, as shown in 
Fig. 9. While there is a complete theory of the statistics 
of dynamics in reflection, in transmi ssion only the ID 
case, = 1, has been solved exactly |l40| |. In this case, 
the transmission delay time is found to be the mean of 
reflection delay times of the both ends of the wave guide . 
For > 1, an approximate solution was found [l40j, 
which is shown in Fig. 10a, and compared to results of nu- 
merical simulations. Microwave measurements of P{4>') 
within the localization window were carried ou t in e nsem- 
bles of alumina samples of increasing length |l34j . The 
distributions of the normalized single-channel delay time 
in samples with L — 2b and 90 cm are shown in Fig. 10b. 
The values of var(sa) in these samples are 1.0 and 7.1, 
and the ensemble-averaged values of (/)' are (0') = 21 and 
122 ns, respectively. As Fig. 10b suggests, P{(j)') becomes 
more asymmetric as L increases, and reaches its peak at 
values of 4>' further below its average value of unity. The 
measured distributions of Fig. 10b appear to be similar to 
the calculated distributions in Fig. 10a, though a closer 
comparison reveals a discrepancy in the tail that may be 
due to the presence of absorption in alumina samples. 

The averaged delay time is found by summing the 
single-channel delay time weighted by the intensity, 
Wab = lab'P'ab- When Summed over all input and out- 
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FIG. 9: Calculations of distribution of the single-channel re- 
flection delay time in th e dif fusion (top panel) and localized 
(bottom panel) regimes |l32j| . The maximal value of -P (</>') 
in the case a = a' of equal excitation and reflection modes 
(open circles) is greater than that for a ^ a (solid circles) by 
a factor close to \/2. The effect appears only in the localized 
regime, (with kind permission of the authors) 

put modes, Wab relates directly to a fundamental dy- 
namic quantity in condense matter, namely the number 
of microstates per unit fre quen cy interval du inside the 
scattering medium, N{lli), |141| 

-KblabC^'ab = N{CJ). (5) 
TT 

The summation runs over N channels in both reflec- 
tion and transmission. Though the density of states 
in an open (scattering) system is ill-defined, the l.h.s. 
of the Eq. (5) is wel l defined and proportional to 
/ c?r — l). I141| which equals the stored electro- 

magnetic energy |l42j . In the diffusive limit of Gaussian 
field statistics, th e distribut ion of = Wab is a double- 
sided exponential |l3nHl3 ij . In a ID sample with L > ^, 
an algebraic decay with exponent -4/3 was calculated for 
P{W) using a model of resonant transport via localized 
states 136] . The probability distribution of the normal- 
ized weighted delay time, p(w= W/(W)^ , for local- 
ized waves in a quasi- I D alu mina samples with L = 90 
cm is shown in Fig. 11 |l34j . The distribution P{W) is 
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FIG. 10; (a) Distribution of the single-channel transmission 
delay time obtained from numerical simulations of random 
scattering in planar waveguides with N = 2 (open circles) and 
N = 30 (filled circles) channels . Th e solid line is the analytical 
prediction, independent of A'' |140|I . (with kind permission of 
the author); (b) Probability distribution of the normalized 
delay time, P{4''), in alum ina samples with L — 25 (circles) 
and 90 cm (squares) Il3 



extraordinary broad, reflecting the enhanced probability 
of long dwell times and the increased variance of dwell 
times at larger values of the intensity for localized waves. 
It can be well approximated by a double-sided stretched 
exponential to the power of 1/3 (Fig. 11). The affect of 
absorption on the form of the measured distributions was 
not studied. 

Notice that, though less probable, negative "delay" 
times for the multimode case are observ ed and are also 
allowed by scattering theory |l43l Il44| . Negativ e de - 
lay times arise as a result of pulse reshaping |l24l| . 
This is in contrast to the "proper" delay times defined 
as the eigenvalues of the Wigner-Smith matrix Q — 
-iS* ■ dS/duj and its trace Tr Q = Y.ab ^ab = 7rN{uj) = 
th, called the Heisenberg time th] the channel- average 
th/2N = Tw is associated with the Wigner-Smith de- 
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FIG. 11: Probability distribution of the normalized weighted 
delay time, P{W), for localized waves (circles) in alumina 
sample with L = 90 cm. The solid line is the model distribu- 
tion, P{W) = aexp {-b\W\^^'^), with a = 0.44 and b = 2.42 

for W > 0, and a = 0.07 and b = 5.50 for W < 0. From 
Ref. Il3l. 



lay time tw . Substantial theoretical effort in mesoscopic 
physics has focused on the probability distributions of the 
"proper" delay tim es and the Wigner-Smith delay time 
|l45lll46lll47L|]'4^ . whose measurements require knowl- 
edge of the entire Wigner-Smith matrix. The joint dis- 
tribution of the inverse "proper" delay times, 7„ = 1/t„, 
was fo und t o be the Laguerre ensemble of random-matrix 
theory |l4(ij , in which the correlation functions of r„ 's are 
given in terms of Laguerre polynomials, whereas P(tw) 
was found to follow univers al distribution s with alge- 
braic decays in the localized 

[m El 1113 

and critical 

|l49l llSOjl regimes, respectively. It has also been shown 
that the statistics of delay times relates to the statistics 
of the eigenfunction intensities and may be used to probe 
eigenfunction fluctuation s in close d mesoscopic samples 
of any spatial dimension |l5ll Il52| . 



STATISTICS OF PULSED TRANSMISSION 

In addition to probing the time evolution of aver- 
age transmission and its increasing suppression, and the 
statistics of delay time for incident pulses, it is possible 
to measure the spatial and polarization correlation func- 
tions and the f ull di stribution of the fleld and intensity at 
any value of t jl53l |. The cumulant correlation function 
with displacement and polarization rotation of intensity 
of the transmitted wave normalized to its ensemble aver- 
age value for an incident pulse of linewidth a at time t 
has the form, 



C = F + KAt)(l + F) 



(6) 



13 



where F is the complex square of the field correlation 
function. The parameter Ka{t) establishes the degree of 
correlation at a delay time t from the center of an excit- 
ing pulse of bandwidth cr, and is evidently the fractional 
intensity correlation at a point at which F vanishes. This 
is the same form as the corresponding steady-state corre- 
lation functions with th e degree o f steady state correla- 
tion kq instead of Ka{t) |l54 llSiil l. The field correlation 
function is found not to be a function of time, while the 
degree of intensity correlation (t) varie s with time and 
depends upon the pulse bandwidth |l53j |. This is seen in 
the polarization measurements shown in Fig. 12 carried 
out in an ensemble of 12,000 configurations of low den- 
sity alumina spheres at L = 90 cm over the frequency 
range 16.95-17.05 GHz near the peak of t he fo urth Mie 
resonance with spectral width a — 5 MHz |l53j . For the 
conditions given above, kq = 0.29. The dynamic proba- 
bility distribution of normalized intensity is found to de- 
pend only upon K„{t), which is equal to the variance of 
the normalized total transmission at a fixed delay time, 
vaT{Sa{t)). The intensity distribution at a fixed delay 
time is given by Eq. (1) with g = 2/3k and closely cor- 
responds to the measured intensity distribution, even for 
K ^ 1. Values of Ka{t) much greater than kq are found 
at long delay times. For an incident pulse width a — 1 
MHz in a sample with L = 61 cm, in the frequency range 
9.95-10.05 GHz, just above the first Mie resonance, at 
t = 740 ns, we find k = 9.55. 



FIELD STATISTICS 

The isotropic approximation of random matrix theory 
is equivalent to the assumption that the field distribu- 
tion in a given configuration is a circular Gaussian |l56| . 
This is demonstrated by measurements of the microwave 
field in random media, which show that the probability 
distribution, P{E), of the field normalized to the ensem- 
ble average of the field amplitude, E — tab/\/ (|ia6p) = 
r+\/ —li, is given by substituting 



r'^+i'^ for Sab in Eq. (1) 



|l57i| . The probability distribution for the transmitted 
field for a sub-ensemble with a fixed value of Sq is thus 
Gaussian. Because the field in samples with a fixed value 
of Sa is a Gaussian random process, the field normalized 
to its average value in a given configuration and the to- 
tal transmission are statistically independent. The field 
correlation function may therefore be written as a prod- 
uct of correlators of the normalized Gaussian field and 
the square root of the total transmission. Since the field 
correlation function with frequency shift is the Fourier 
transform of the timc-of- flight distribution [2^ l30| , pulse 
dynamics, which shows progressive suppression of trans- 
port, may be expressed in terms of the spectral correla- 
tion function of the normalized Gaussian field and the 
square root of total transmission. We expect that the 
correlation function of the Gaussian field will not be af- 
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FIG. 12: (a) Real part of the field correlation function, 
Fe{A9), and (b) intensity correlation function, Ccr{A9), with 
polarization rotation of the transmitted wave at the two de- 
lay times following pulsed excitation by a Gaussian pulse with 
spectral width cr = 5 MHz and for monochromatic excitation 
(GW). The dashed curves are (a) FsiAe) = cos(A6') and (b) 
Ca = F{Ae) + Ka{t) [F{Ae) + 1], where F(A6») = \FE{Ae)\^ 
with the values of Ko-(t) indicated. The solid curve is F{A6). 
The logarithm of the average pulsed transmission through the 
alumina sample for a = 5 MHz, normalized by the average 
steady-state transmitted intensity, is shown in the inset. The 
delay times at which correlation is presented i n the figure are 
indicated by vertical dashed lines. From Ref. |153| . 



fected by mesoscopic fluctuations. On the other hand, 
the correlation function of the square root of the total 
transmission reflects mesoscopic fluctuations. In this way 
the dynamics of the field with the observed suppression 
of transmission in the time domain due to the growing 
impact of localization is directly tied to mesoscopic fluc- 
tuations. 



CONCLUSION 

Because the character of wave transport depends upon 
the proximity to the localization threshold, signatures of 
localization, and of the approach to localization, can be 
found whenever statistical measurements can be made in 
ensembles of random sample realizations. A large num- 
ber of statistically different sample configurations may 
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arise as a result of constant motion of colloidal systems, 
or in powdered samples, or in collections of scattering 
particles that can be tumbled. They may also be cre- 
ated by illuminating different parts of a large slab or by 
tuning the frequency of a monochromatic source through 
a spectrum much wider than the correlation frequency. 
The probability distributions of intensity and total trans- 
mission in steady-state and pulsed measurements can be 
described in terms of the variance of normalized total 
transmission, which is equal to the degree of intensity 
correlation, var (sa) = k. The variation of k with delay 
time directly charts the increasing impact of localization 
and is independent of absorption. The decrease of fluc- 
tuations and correlation with absorption in steady-state 
measurements reflects the reduced weight with increas- 
ing absorption of long trajectories associated with strong 
correlation. 

Absorption reduces both steady-state and pulsed 

transmission. The influences of absorption and weak lo- 
calization can be isolated in the time domain. Absorp- 
tion reduces the intensity of pulsed transmission expo- 
nentially. This reduction can be cancelled in the analysis 
of data by multiplying by the corresponding growing ex- 
ponential, to yield the transmitted pulse that would be 
observed in the absence of absorption. The role of lo- 
calization in averaged transmission at a given time delay 
is thus displayed directly. In contrast, absorption and 
localization affects are mixed in steady-state transmis- 
sion, and either can lead to an exponential decrease of 
transmission with sample thickness. 

The suppression of the spread of a beam within a sam- 
ple also gives an indication of localization. This can be 
observed in the broadening of the peak and in the reduc- 
tion of the coherent backscattcring enhancement factor 
below 2. Other signatures of localization are observed in 
measurements of the statistics of the single channel de- 
lay time and its correlation with the intensity. Given the 
many signatures of localization, the state of propagation 
can be unambiguously determined even in the presence 
of absorption. 
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